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Effects of inhomogeneous doping on the high-T c cuprate superconductors are studied within the 
framework of the t-J model. Especially, the boundary between two non-superconducting regions 
with doping rates much higher and lower than the optimal one is examined. It is found that, although 
there is no superconductivity in the bulk, a superconducting region appears at the boundary because 
singlet resonating valence bond order and holon condensation occur simultaneously in this region. 
The critical temperature of the induced superconductivity can be higher than that of the optimally 
doped sample and the critical current density can be the same order as that in the superconducting 
state below the critical temperature. We also point out an experimental possibility to observe this 
phenomenon. 
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It has been recognized that there are two important in- 
gredients for the superconductivity in high-T c cuprates, 
which are the spin correlation, leading to the spin gap, 
and the doped charge carriers with a sufficiently high 
density. In the slave-boson mean- field theory of the t-J 
model, the former is realized as the singlet resonating va- 
lence bond (s-RVB) order of spinons and the latter the 
bose condensation of holons (l[ ||, ||, |[ [|, Q . In this 
framework, it is assumed that the spinons and the holons 
are well-defined low lying excitations, and the spin-charge 
separation is a good picture. Although there still remains 
some controversy on this point, it has been clarified that 
some physical properties of the high-T c cuprates can be 
well described by the t-J model, for example, the trans- 
port properties || and the magnetic properties (|, |To| . 
In this Letter, we take the t-J model to be our starting 
point. 

We consider the situation where the doping rate is spa- 
tially varying. We especially concentrate on the case 
where two non-superconducting regions, one is too over- 
doped and the other is too under-doped, are in contact 
with each other within the same Cu02 plane. We expect 
that the s-RVB order and the holon condensate can ex- 
ist simultaneously at the boundary, thus giving rise to 
superconductivity. Since the s-RVB order and the holon 
condensation can occur at temperatures higher than the 
critical temperature of the optimally doped sample, T opt , 
this "boundary superconductivity" may also occur above 

Our idea is similar with the "spin-gap proximity 
effect" previously introduced by Emery, Kivelson and 
Zacher based on a somewhat different model for high 
T c superconductors, discussing the stripe order and its 
significance for the superconductivity. Recently, Kivel- 
son has proposed a possibility to raise critical temper- 
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ature using spin-gap proximity effect p2[ . He proposes 
making a stack of spin-gap material and hole-rich mate- 
rial in a layered structure. In contrast to this, we consider 
the situation where the spin-gap and hole-rich regions ex- 
ist within a Cu02 plane. This kind of situation has not 
been studied within the framework of the t-J model and 
we consider this is of interest from both theoretical and 
experimental points of view. 

First we introduce the t-J model which describes the 
correlated electrons in each Cu02 plane. The Hamilto- 
nian is given by 

« = -*E (flfiabfa + h-c) + J E §i ■ Sj 

-E A *(E^-+ 6 ^- 1 )' w 

where i and j denote lattice points in the Cu02 plane, 
means the nearest neighbors, fi a (a =j, j) and 
bi are the annihilation operators of the spinon and the 
holon, respectively, Si denotes - ^ a(3 fj a a af3 f if3 with a af3 
being Pauli matrices, and A^ is the Lagrange multiplier 
for the constraint, ^2 a fi a fia + b\bi = 1. In this Letter 
the lattice spacing is taken to be unity. In addition to 
Eq. (|l|), we later introduce the coupling to the electro- 
magnetic field. 

The Eq. (|l|) is treated by the so-called slave-boson 
mean-field theory. The local constraint represented by 
A^ is replaced by the global one, which is described 
by the constant chemical potential for spinons /ip and 
that for holons ji^. We define the statistical averages, 
£ = (J2 a fiafja) an d tj = (b]bj), which are assumed to 
be independent of i and j Q . The order parameter for 
the s-RVB order is denoted as A^- = jfiftffji — Ai/jt)? 
which is equal to A x (fi) ifj = i+x and A y (fi) if j = i+y 
where x and y indicate the shift by one lattice spacing in 
positive x- and ^-direction, respectively. Although £, r] 
and Aij should be determined self-consistently in princi- 
ple, we skip this process for simplicity by just assuming £ 
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and 77 to be constant. Including also the scalar and vector 
potential, ip and A, and the fictitious gauge fields, ao and 
a, which represent the fluctuation around the mean-field 
solution ||], we obtain the following Lagrangian, 

£ = ^2fi* i hd r - + iaoin)} fia 
ia 

+ ^2 b i { hdr ~ Mb + i(ecp(fi) + a (fi))} h 

i 
i 

-(^ + ^)E{ e< -/U.+h.c.} 

-«£{e«36l6 J + h.c.} 



A*(ri)xii+£ + A*(ri)xii+y + h.c. 



^(lA^rD^ + IA^rl)! 2 ) 
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where Xij = fi]fjl ~ fnfj] and 



^=U/-^ e *=\!?S a+ -/)- Al (3) 

The chemical potentials /ip and /ie are determined from 

= 1-1 (b\bi) = 5, (4) 

where 5 is the average holon density in the whole system. 
The background charge density, originating from doping, 
is denoted by p(fi), whose spatial average equals — eS. 
In this Letter ao and a are treated perturbatively on the 
same line with Refs. || and Q. 

First we study how the distribution of the holons, 
namely 5(r£), is determined under an inhomogeneous 
doping, which is described by p(fi). For this purpose, 
we neglect the terms including A^ and a and A. Then 
the spinon and the holon degrees of freedom can be in- 
tegrated out, giving rise to terms, ttfclo + ttb^o + e{ PY \ 
where 7Tf and 7Tb , respectively, are the polarization func- 
tions of spinons and holons. Now the magnitudes of 7Tp 
and 7Tb determine which of the spinon and the holon cou- 
ples to cp. In this Letter we assume 7Tp ^> ttb, which may 
be valid in lower doping region. Then ao is almost fixed 
to zero and (p couples to the holons. Now S and cp are 
determined from p in a self-consistent way. 

We especially consider the following situation. The 
Cu02 planes, described by the above Lagrangian, are 
laid parallel to x-y plane and stacked in ^-direction. The 
region x < (x > 0) is filled with a s-RVB ordered 
(hole-rich) material with a low (high) enough doping rate. 
Both of the regions are non-superconducting. The system 
is uniform in y- and z-direction. We assume that the 
doping rate, which is described by p(ri), changes at x = 
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FIG. 1: Spatial variation of the holon density S (dotted line) 
and the s-RVB order parameter (bold line). H and L 
stand for high and low doping region, respectively. In the 
shaded region the superconductivity is expected to appear. 



abruptly. In order to preserve the charge neutrality in 
the bulk, namely at x — > ±00, the statistical average 
(<p(r)) must take constant imaginary values, which differ 
between x — > 00 and x — > — 00. (Note that the chemical 
potential pb is a constant in the whole system.) Near the 
boundary (ip(r)) changes smoothly connecting these two 
bulk limits. As in the case of semiconductor junctions 
]l5[ , the smooth connection is enabled by the polarization 
charge appearing at the boundary. The depth of the 
charged region is given by the screening length of the 
holons. This length is extremely long (short) for x < 
(x > 0) and there still remains a sharp drop of the holon 
density at the boundary. Therefore, in this Letter, we 
approximate 5 by a step function (see Fig. [I]). 

Next we study the s-RVB order by introducing the 
Ginzburg-Landau (GL) free energy of the order parame- 
ter Ad = \(A X — A y ). Here we assume that the s-RVB 
order has d-wave symmetry. The GL free energy is ob- 
tained by the perturbative expansion of the free nergy 
with respect to A^ Jl^| . We introduce the half band 
width D = 4(trj + 3J£/8) and the strength of attrac- 
tive interaction V = 3 J/8. The dispersion of spinons, 
(D/2)(cos 

k x + cosfcy), is approximated by a parabolic 
one, D(—l + |/c| 2 /4), and then the density of states at 
Fermi energy becomes N(0) = 1/(ttD). At a temperature 
T C D, the solution for ^2 a (fi a fia) = 1 — (5 is approxi- 
mately given by p^ = —S/N(0) and the critical temper- 
ature T d becomes T d ~ T* exp[-l/{2WV(0)(l - 5) 2 }} 
with T* = 2De 7 /7r where 7 is the Euler's constant. 
This expression is valid in the "weak coupling regime", 
VN(0) <C 1. The GL free energy density of A^ is given 

by 
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a d \A d \ 2 - 



/3|A d | 4 + 7d |nA d 
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where II = -iV+ (2/H)a+ (2ir/<p )A with <p = hc/(2e 
The GL coefficients are given by 

T 



ad 

P 
Id 



2iY(0)(l-<5) 2 ln 



21C(3)JV(0) 



(1-5) 4 



27T 2 T 2 v > - T 2 D > 
7C(3)7V(0)£> 2 = C2D 
16tt 2 T 2 ~~ T 2 ' 



(6) 
(7) 
(8) 
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where ((x) is the zeta function and c\ and C2 are con- 
stants. The doping dependence of the GL coefficients 
comes from X^, otd and (3. Since the largest dependence 
appears from X^, we neglect the (1 — (5) n -factors in ad 
and (3 in the following. Then the doping rate affects Td 
only. 

Here we consider the situation where the transition 
temperatures in the regions x < and x > 0, denoted 
by Td and T^, respectively, satisfies the condition, Td > 
T ^> T' d . It is also assumed that holons are condensed in 
the region x > 0. Note that the quantity with (without) 
dash (') is defined in x > (x < 0). We first determine 
the spatial variation of A^, disregarding a and A. Here 
Ad is assumed to be real. The expectation value of A^ 
is y/\ad\/(2f3) = Ad at x — > — oo and at x — > oo. In 
x < we introduce a new variable JA^ = A^ — A^. The 
spatial variation of A^ is then governed by the following 
equations, 



2\a d \5A d (x) - <y d %5A d (x) =0 (x < 0), 
a' d A d (x) - ldd 2 x A d {x) =0 (x > 0). 



(9) 



We can easily see that A^ behaves as A^ + gi exp(x/^p) 
for x < and #2 exp(— #/£p) for x > 0, where £f = 
V^7d/2|ad|, = y/jd/a'd, and #i and # 2 are constants. 
The solution in the whole region is obtained by connect- 
ing these two solutions at x = with boundary condition, 

d x A d (-0) = d x A d (+0), Ad(-0) = A d (+0), (10) 

(see Fig. Q). This condition is in contrast with that for 
the proximity effect in ordinary superconductor-normal 
metal junction, where the change of the BCS interaction 
causes a discontinuity of the gap function [[l7|. (Note 
that the quantity which is continuous at the boundary is 
not the gap function but the anomalous amplitude.) The 
constants g\ and #2 are obtained as 
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£f + £ f 
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A d . 



(11) 



Note that non-zero Ad in x > is due to the proximity 
effect of s-RVB order and £ F gives the proximity length. 
The temperature dependence of the coherence lengths at 
X^ <C X <T d are as follows: 



c 2 D 



1 



2 Tdy/l-T/Td 



D 



It is interesting to see that £ F is, except for the logarith- 
mic factor, analogous to the ordinary proximity length in 
the clean limit ex hv^/T where is the Fermi velocity. If 
T is close to T d , £f becomes much larger than £ F and the 
amplitude of Ad in the proximity region, i.e., #2, becomes 
much smaller than A^. However if T is not too close to 
Td, £f and £ F can be the same order and #2 becomes com- 
parable to Ad. Noting that the holons are already bose 
condensed in the proximity region, this result means the 
induced superconductivity in the boundary area. 



The strength of the induced superconductivity be- 
comes more apparent by estimating the supercurrent 
which can flow parallel to the boundary. In this region 
the gauge field a is massive because of the holon con- 
densation, and a in Eq. (j^j can be set to zero. Then the 
response to the actual vector potential A{r) is dominated 
by the spinons. Here we take the vector potential to be 
constant and parallel to ^/-direction. The total current in 
^/-direction is given by 



J 



r>00 

2c ld A y / A 2 d (x)dx 
Jo 



= CJdAyJJJ 
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where the A^-dependent proximity length £p(A y ) — 
£ F / ^Jl + p' A 1 * decreases with increasing A y , with p r be- 
ing (27r£ F /</>o) 2 . The maximum of J as a function of A y 
is given approximately by 

2 



/ 3 




(13) 



For clarity, here we introduce the critical current density 
j c in the bulk superconducting phase. This quantity (in 
overdoped region) is calculated from Eq. (|J) by assuming 
that holons are condensed and a is massive in the whole 
region. We obtain 



/2tt\ 2 1 ■ 

Jc = 2c {To) 7d vT 



(14) 



where p = 2(27t£f/0o) 2 - It is clear that J ma x = 
V^JcCf 2 /^- If £f — £f 5 the critical current density in the 
proximity region, defined by Jmax/^F? can be the same 
order as j c . 

Finally we point out some experimental set-ups to ob- 
serve the effect predicted in this Letter. The simplest way 
may be to utilize the field effect transistor (FET) recently 
developed by Shon et al. []l8|, |l9|| . We propose two types 
of experiments, which we call "half-gate" experiment and 
"side-gate" experiment. The experimental set-ups are de- 
picted in Fig. H In the former case (Fig. ||(a)) the gate 
electrode is attached parallel to the Cu02 plane, but it 
covers only half of it. As a result, a boundary between 
the s-RVB ordered region (uncovered by the gate) and 
the hole-rich region (covered by the gate) is formed un- 
der the edge of the gate. (Here we are assuming that the 
gate induces doping to the undoped compound.) How- 
ever, in order to observe the effect studied in this Letter, 
the insulating layer separating the gate and the sample 
must be smaller than the s-RVB coherence length (ap- 
proximately the superconducting coherence length below 
T c ~ several nm.), which may be hardly realized. 

The other configuration is the "side-gate" as depicted 
m Fig. |(b). In this case the gate is applied to the Cu02 
plane from the side and by applying a voltage a hole- 
rich layer is formed under the gate. This situation is a 
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FIG. 2: Experimental set-ups using field effect transistor: (a) 
half-gate configuration and (b) side-gate configuration. The 
Cu02 planes are shown in an exaggerated scale as compared 
to the electrodes. 



little different from one considered in this Letter where 
the hole-rich region is infinitely wide. However the same 
mechanism studied by us may also work. In this case, the 
thickness of the superconducting region is limited not by 
the s-RVB proximity length but by the thickness of the 
hole-rich region. This configuration can be realized more 



easily. 

There may be other ways. One might introduce the 
spatial variation of doping chemically, for example, by 
using the epitaxial growth |20| . If this can be realized, it 
may fit the situation studied in this Letter best. 

In either case, the width of the superconducting re- 
gion is very narrow and the fluctuation is essential 
Therefore it is not easy to observe good superconducting 
behavior. However we expect that whether the super- 
conducting correlation exists or not is an experimentally 
detectable fact. 

In summary, we have studied the effect of inhomoge- 
neous doping in the high-T c superconductors. We have 
shown, based on the t-J model, a possibility of supercon- 
ductivity at the boundary between the s-RVB ordered 
and the hole-rich regions, caused by the proximity effect 
of the s-RVB order. We have also suggested experimen- 
tal set-ups to observe this effect. Since our proposal is 
based on the spin-charge separation picture, these exper- 
iments may also shed a renewed light on the mechanism 
of high-T c superconductivity. On the other hand, from 
a theoretical point of view, the extension of the present 
study to more general inhomogeneous cases is also im- 
portant in understandin g r ecently found inhomogeneous 
phenomena in cuprates |22], ^3|, |24| . 
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